A hyperboloidal study of tail decay rates for scalar and Yang-Mills fields 
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We investigate the asymptotic behavior of spherically symmetric solutions to scalar wave and 
Yang-Mills equations on a Schwarzschild background. The studies demonstrate the astrophysical 
relevance of null infinity in predicting radiation signals for gravitational wave detectors and show 
how test fields on unbounded domains in black hole spacetimes can be simulated conveniently by 
numerically solving hyperboloidal initial value problems. 



I. INTRODUCTION 

The purpose of this article is two-fold. First, we em- 
phasize the astrophysical relevance of null infinity in pre- 
dicting gravitational radiation signals as they are ex- 
pected to be detected by gravitational wave detectors. 
It has been demonstrated in [40] that a suitable feature 
to estimate the validity of the concept of null infinity 
when predicting detector signals is the tail behavior. In 
this article, we work out this idea further by discussing 
the observer dependence of tails for solutions to scalar 
wave and Yang-Mills equations on a Schwarzschild back- 
ground. Secondly, it is demonstrated that test fields at 
null infinity can be conveniently studied numerically by 
solving Cauchy problems once a suitable choice of gauge 
for the background spacetime has been made. 

The question about the astrophysical relevance of null 
infinity is a conceptual one. Mathematically, it is only at 
null infinity that one can give an unambiguous definition 
of gravitational radiation. Many physicists, however, are 
uncomfortable with the notion of null infinity because the 
detectors are at a finite distance away from the source 
and move along timelike curves while idealized observers 
at null infinity are infinitely far away from the source 
and move along null curves. These issues are addressed 
within the article on concrete examples. 

Our aim is not to discuss the gauge dependency of radi- 
ation extraction or the influence of finite distance bound- 
ary conditions on the numerical solution as has been done 
for example in (35l . l4l| . We will rather numerically study 
a genuinely physical phenomenon, namely the crucial de- 
pendence of a certain feature in the radiation signal on 
the observers distance from the source [24j, |40( • 

There is compelling evidence suggesting that, asymp- 
totically with respect to an appropriate time coordinate, 
the gravitational radiation flux from a p erturbed black 
hole decays polynomially in time [H, OjJ [H, HE [2(| [3^| . 
This polynomial decay can be regarded as being due 
to back-scatter of gravitational radiation off the curved 
background spacetime. The back-scattering effect is very 
small but it may, in principle, be detected by future gen- 
eration gravitational wave detectors [f| . The decay rate 



'Electronic address: anil@aei.mpg.de 



of the tail may be of particular interest in testing funda- 
mental predictions of general relativity because unique- 
ness theorems suggest that all features of black holes 
other than mass, charge and angular momentum will de- 
cay according to this rate independent of the details of 
the collapse process. The tail plays an important role also 
in the mathematical discussion of stability of black hole 
spacetimes and the cosmic censorship conjecture (l5lll6j. 

It turns out that the gravitational radiation flux of 
a general perturbation has a faster decay on the event 
horizon and on timelike surfaces than on null infinity 
[1 H 0, 03, E HE HI. This difference raises the ques- 
tion of which decay is the relevant one with respect to our 
gravitational wave detectors. Piirrer et al have demon- 
strated in [Z(| that the relevant decay rate in the so-called 
astrophysical zone [HJ is the one measured at null infin- 
ity. Their study is based on a characteristic approach 
and discusses critical collapse for self-gravitating mass- 
less scalar fields in spherical symmetry. 

We are interested in a detailed study of the observer 
dependence of decay rates including null infinity and as- 
sume the background to be given. We compare a pow- 
erful mathematical estimate on decay rates for scalar 
fields suggested in [HI, [43[ with numerical calculations 
and check decay rates for scalar fields with non- vanishing 
angular momentum. We also study the tail behavior for 
Yang-Mills fields where a breakdown of linear perturba- 
tion theory has been demonstrated at future timelike in- 
finity 5] . Our study indicates a similar breakdown of the 
linear analysis along future null infinity. 

Additionally, we show how test fields on black hole 
spacetimes can be studied numerically by including null 
infinity while solving Cauchy problems. Accurate studies 
of fields on black hole spacetimes including null infinity 
have been restricted to the characteristic approach until 
now [l2|, HH, (HI- The characteristic approach, however, 
has certain difficulties mainly due to problems related to 
the rigidity in the choice of the underlying gauge [Zij . 
The method presented in this article is quite flexible and 
allows us to use standard discretization techniques com- 
monly applied in numerical relativity. The feature that 
makes the accuracy of our calculations possible is that 
the numerical outer boundary coincides with null infinity 
so that, in the continuum problem, no boundary condi- 
tions are required and no resolution loss appears in the 
physical part of the conformal extension [ll 1X71. l30l. l4r| . 
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II. THE CONFORMAL METHOD 

We want to study solutions to the scalar wave equation 
□ $ = 0, (1) 

for a scalar field $ in a Schwarzschild spacetime. To in- 
clude null infinity in the computational domain in a man- 
ifestly regular fashion, we apply the conformal method in- 
troduced by Penrose [a^Halll (see also [H IM HI HI ) • 
A conformal rescaling g = fl 2 g of the physical metric g 
with a conformal factor Q > implies the transformation 



□ - -R ) $ = fT 3 
6 



1 ~\ - $ 
□ --i?U, with $=-, 



where R and R are the Ricci scalars of the rescaled and 
the physical metrics g and g respectively. The equation 
(JXJ) with respect to the rescaled metric then becomes 



□ $ - -R<f> = 0. 
6 



(2) 



We are also interested in solutions to the Yang-Mills 
equations, 



D*F :=d*F + Af\F = 0, 



(3) 



where F = dA + A A A is the Yang-Mills curvature and A 
is the Yang-Mills connection. The Yang-Mills equations 
are conformally invariant. We can therefore study the 
system $5$ directly in a conformally rescaled spacetime. 

Schwarzschild spacetime is weakly asymptotically sim- 
ple, implying that a suitably rescaled metric is smoothly 
extendable through null infinity. We can solve the equa- 
tions and ([3]) in a conformally extended Schwarzschild 
spacetime by allowing the conformal factor to vanish in a 
certain way provided that the initial data has sufficiently 
fast fall-off [23. 

To solve an initial value problem for ([2]) and ([3]) numer- 
ically on a Schwarzschild spacetime, we need to choose a 
coordinate system that is suitable for numerical calcula- 
tions and covers that part of the Schwarzschild-Kruskal 
manifold we are interested in, namely the domain extend- 
ing from a neighborhood of the event horizon up to future 
null infinity. Our coordinates should also be adapted to a 
natural family of observers, each being thought as repre- 
senting gravitational wave detectors at constant distances 
away from the black hole. 

A convenient foliation covering the extended Schwarz- 
schild spacetime up to future null infinity satisfying the 
above requirements can be found in the class of spher- 
ically symmetric constant mean curvature (CMC) folia- 
tions of Schwarzschild spacetime @, H2, [HI . The trans- 
formation from the standard Schwarzschild time coordi- 
nate to the time coordinate of a spherically symmetric 
CMC-foliation can be written as 



t = t-h(f), 



(4) 



where f is the standard Schwarzschild area radius, h(f) 
is the height function and t is the standard Schwarzschild 
time coordinate. The derivative of the height function is 
given by 



h'(f) 
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C 



(1 



p(fy 



(5) 



with 



P(f) := 



\ 



Af 3 



-C 



2m 
f 



The mass of the Schwarzschild black hole is denoted by 
to. The foliation parameters are the mean extrinsic cur- 
vature A, and a constant of integration C . The global 
behavior of CMC-surfaces depends on the foliation pa- 
rameters. We choose the parameters such that the sur- 
faces of the foliation come from future null infinity, pass 
the event horizon above the bifurcation sphere and run 
into the future singularity as depicted in Fig. [T] 

singularity 




Fig. 1: Penrose diagram of a CMC-foliation in Schwarzschild 
spacetime with m = 1/2, C = 1, K = 2. The dashed lines 
represent Killing observers. 

The transformation (J4j) implies that the coordinates of 
the foliation are adapted to time symmetry in the sense 
that the timelike Killing vector field is given by dt- To 
keep this property for the conformally rescaled metric, 
we perform the conformal compactification such that the 
conformal factor is time-independent which we refer to as 
conformal fixing. The Schwarzschild metric in a CMC- 
foliation and a conformal fixing gauge can be written as 

m 



g = - I 1 - — ) 



■ dr 2 + r 2 da 2 



(A> 3 /3 



P(r) 



P 2 {r) 

where da 2 is the standard metric on the unit sphere 



dtdr + 
(6) 



P(r) 



\ 



Kr 3 



cn 3 



1 



2m(l 



n 2 r 4 
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and fi = 1 — r. Note that conformal fixing on a hyper- 
boloidal foliation implies scri-fixing, that is, the coordi- 
nate location of null infinity in a conformal fixing gauge 
is independent of the time coordinate. The compactify- 
ing coordinate r is related to the standard Schwarzschild 
area-radius f via 



r 



(7) 



The timelike curves with constant spatial coordinates 
can be regarded as worldlines of natural observers at con- 
stant distances away from the source. As can be seen 
from Fig. Q] and from g(dt, dt)\ y+ = 0, the worldline of 
such observers becomes null at null infinity. 



III. THE INITIAL VALUE PROBLEM 
A. The continuum problem 

We write the metric © as 

g= (-a 2 + 1 2 2 ) dt 2 + 2 7 2 (3 dt dr + 7 2 dr 2 + r 2 da 2 . 

The lapse a, the shift (3, and the spatial metric function 
7 read 




We apply the method of separation of variables to 
We write $(t, r, ti, if) = 0(t, r)Y lm (i?, if) with Y lm being 
the usual spherical harmonics and define the auxiliary 
variables 

V>:=dr0, and tt := ~(d t d> - 0d r d>). (9) 
a 

Then the following linear, symmetric hyperbolic system 
of evolution equations for (f>(t, r) is obtained 

a 

d t 4> = — n + pip, 
7 



is given by = 0. The auxiliary variables are defined 
as in ([9]). The only difference to the system (fTOj) is the 
equation for the time derivative of 7r which becomes 



d t 1p = d r I — 7T 

\7 

1 



(10) 



The Ricci scalar of the metric ([6]) reads 
12ft 



" " ( \ R( t>- 



1(1 + 1) 



R 



(r + m(2r - 1)) 



For the Yang-Mills equations we consider the gauge 
group SU (2) and make the following spherically symmet- 
ric ansatz for the connection 0, fl4j | 



A 



+ 1) a% d-d + (cos ■& a 3 + (<j> + 1) sin a 2 ) d<p, 



where the Ui are Pauli matrices and <f> = <p(t,r). With 
the above choice, a vacuum state for the Yang- Mills field 



dfTT = d r 



-V + /37T) -^(1 + ^(2 + 0). (11) 

7 / r z 



We choose an approximately outgoing, compactly sup- 
ported Gaussian pulse for initial data. On the initial 
hypersurface such data satisfies 

dt[r<j>) +c+dr{r^) =0, 

where c+ is the outgoing characteristic speed given by 
c + = — p + a/j. The data reads for r < r out = const. < 1 
as follows 



0(0, r) 
V>(0,r) 

7r(0,r) 



a e 



2 (r - r c ) 
a 2 

-V(0,r)- 



0(0, r), 
0(0, r) 



r \ a 



(12) 



where r c is the center of the Gaussian pulse, a is its am- 
plitude and er is its width. The data is set to zero for 
r > r out . 

The initial hypersurface is hyperboloidal, meaning that 
it is a smooth spacelike hypersurface extending through 
null infinity. The problem consisting of initial data given 
on a hyperboloidal surface and evolution equations is 
called a hyperboloidal initial value problem [20]. Note 
that no boundary conditions are needed in the contin- 
uum problem. 



B. The numerical method 

We may solve hyperboloidal initial value problems for 
scalar wave and Yang-Mills equations using standard nu- 
merical techniques. I used the method of lines with fourth 
order Runge-Kutta time integration and finite differenc- 
ing with sixth order accurate stencils. The inner bound- 
ary is chosen to be a spacelike surface inside the event 
horizon so that one can use an excision technique. The 
outer boundary is at null infinity where one-sided finite 
differencing is applied. The lack of outer boundary condi- 
tions makes the numerical implementation simpler than 
in the case where an artificial timelike outer boundary is 
introduced. 

The convergence of the code can be seen in Fig. [2j For 
this plot, a three level convergence analysis has been per- 
formed with 1010, 2020 and 4040 grid cells. The Courant 
factor At/ Ar is chosen to be 0.5. The convergence in the 
L2-norm has been plotted for a solution to the hyper- 
boloidal initial value problem for the scalar wave equa- 
tion. The corresponding plot for the Yang-Mills system 
is similar. We see that the code switches from fourth to 
sixth order convergence after a short time. This seems to 
be due to the initially dominating error from the fourth 
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Q 



Fig. 2: Convergence in the L2-norm for numerical solu- 
tions to (|10[) . The convergence factor Q is calculated by 
n — l ||0'° m -0" led || 

V — 10 S2 ||^me(i_^h«BA || • 



order time integration, and indeed, choosing a Courant- 
factor 10 times smaller results in the dashed curve which 
indicates sixth order convergence from the start. As we 
are interested in long time evolutions, the larger Courant 
factor is chosen. The convergence factors in the maxi- 
mum and the Ll-norms behave similarly. 

For the results presented below the number of grid cells 
is chosen to be 10020. The spatial simulation domain in 
grid coordinates is given by r E [0.499, 1]. The Courant 
factor is 0.2. The parameters of the spacetime and the 
foliation are m = 0.5, K = 0.5, C = 1. 



problem gives 



(C 2 +t + r*)(C 3 +t-n) 



(13) 



where (Ci, C2, C3) are constants, f* := f+2mlog(f — 2m) 
and t is the standard Schwarzschild time coordinate. 
This estimate captures the asymptotic behavior of the so- 
lution near null infinity as well as at finite distances from 
the black hole. The corresponding decay rates can easily 
be read off from (fT3"|) . We want to compare the mathe- 
matical estimate with the numerical calculation with re- 
spect to the prediction of the decay rate. We calculate 
the function 



p(t, r) 



dint 



(14) 



which becomes asymptotically the exponent of the poly- 
nomial decay of the solution. For the comparison, we 
need to write down the mathematical estimate (|13j) in 
the coordinates of our numerical calculation. The coordi- 
nate transformation is given by ((4]). The height function 
in this transformation is not given in explicit form, but 
as we are mainly interested in the asymptotic behavior 
of the solution in the far-field zone we can use a Taylor 
expansion of ([5|) as f — > 00. We get 



h(f) 



2K 2 



4m' 



4m J 



IV. RESULTS 
A. Tails for the scalar wave equation 

The late-time decay of a solution to the scalar wave 
equation is expected to have the form 

lim <f>(r,t) = Ct p , 

t — >oc 

where C depends on r only and p is the decay rate with 
p < 0. For a solution with vanishing angular momentum, 
that is, I — 0, we expect along the event horizon and 
along timelike surfaces a decay rate p = —3 [3!|. At 
null infinity, however, the decay rate is expected to be 
p = — 2 0, [25[. We would like to know which decay 
rate is the relevant one in the context of astrophysical 
predictions and have a quantitative understanding of the 
dependence of the decay rate on the observers location 

m. 



1. Comparison between mathematical and numerical results 

Recently, an optimal pointwise decay estimate has 
been suggested, but not yet proven, for small spherically 
symmetric solutions to nonlinear wave equations with a 
potential [H, |43| • Applying this estimate to our linear 
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Fig. 3: Comparison between the mathematical pointwise esti- 
mate (|13|) and a numerical calculation. Each curve depicts the 
function p(t, R) calculated as in (|14[) for the following values 
R/m = (00, 40000, 20000, 10000, 4000, 2000, 1000, 400, 200, 100). 



The comparison has been plotted in Fig. [3l Each curve 
in the figure gives the function p(t, R) for a constant value 
of R. The dashed lines represent the mathematical esti- 
mate which have been fitted at late times to some repre- 
sentative data points from the numerical simulation. We 
set C\ = C% = 1 and treat C3 as a fitting parameter along 
each curve. Considering that the formula (fl~3|) is only an 
estimate and not an approximation to the solution, we 
can say that it is remarkable how well the qualitative 
behavior of the asymptotic decay is captured by it. 
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2. Astrophysical relevance of null infinity for the decay rate 

Now we can discuss the question which decay rate is 
the relevant one for gravitational wave detectors. We 
think of each curve in Fig. [3] as representing the mea- 
surement made by a detector at some constant distance 
away from the source. In a Schwarzschild spacetime we 
have a length scale at our disposal given by m, the mass 
of the black hole, so that we can gain an intuition for the 
scales by considering the distances of some typical astro- 
nomical objects in terms of geometric units. The closest 
candidate for a supermassive black hole is Sgr A* at the 
center of our galaxy. Its mass is about 3.7 x 10 6 Mq and it 
is about 26000 light years away. In geometric units, this 
corresponds to roughly 10 10 m. Another close example 
is the first serious candidate for a stellar size black hole, 
namely the compact object in the binary system Cygnus 
X-1. The total mass of the binary system is about 50Mq 
and its distance to the Earth is about 8000 light years. 
This gives in geometric units a distance of roughly 10 15 to. 
Now consider in Fig. [3] the curve that is closest to what 
has been calculated along J r+ . Its distance to the black 
hole in geometric units is 4 x 10 4 m. While this is quite far 
for numerical calculations, it is obviously still very close 
compared to astronomical distances. Even for such com- 
parably small distances we see that the relevant decay 
rate is the one measured at null infinity. 

One can of course argue that the decay rate along any 
timelike surface at a finite distance will eventually ap- 
proach —3 if one waits long enough. It is clear from Fig. [3] 
and from the estimate (fT5]> that the time scale for this 
to happen depends on the distance of the source. If we 
have a source that is, for example, 1000 light years away 
from the Earth, then the decay rate will come close to —3 
only after about 10000 years. This is in accordance with 
the argument that the mathematical idealization of the 
region r » t for radiative properties of fields corresponds 
to null infinity while for t ^> f* the correct idealization is 
timelike infinity [TJ, [H [1| . Fig. [3] depicts the transition 
zone between these two regions. We should mention that 
it is extremely unlikely that the polynomially weakening 
tail signal can be followed such a long time that the de- 
cay rate at timelike infinity will be relevant. It is rather 
likely that if ever the tail decay rate is measured by di- 
rect observation, the measurement will be possible only 
for a short time. 

The above discussion recapitulates that is the ap- 
propriate idealization for an observer's location in an 
asymptotically flat spacetime. One should, however, be 
aware that future generation detectors will observe vari- 
ations in the radiation signal due to cosmological effects 
when the sources are considerably far away in compari- 
son to the size of the known universe. In this case, one 
might reconsider the relationship between the different 
idealizations of isolated systems and cosmological space- 
times. 



3. Construction of a physically motivated time coordinate 

In our discussions and conclusions we use the coordi- 
nate time t. In this subsection we want to briefly elabo- 
rate that the same conclusions can be drawn with respect 
a physically motivated time coordinate. 

We have chosen our foliation suitably such that 
worldlines of natural observers of gravitational radia- 
tion are represented by timelike coordinate lines. The 
parametrization of these curves by t, however, can be 
regarded as being arbitrary. There are two steps to be 
taken in the construction of a more physically motivated 
time coordinate. First, we can synchronize our observers 
by an outgoing light pulse so that they have a common 
starting time for their measurements. This starting time 
tn for an observer located at f = R can be calculated by 

pT i-oc 

/ dt = T — tji = / c + (f) dr, 
Jt R Jr. 

where T is the time at infinity and c+ is the outgoing 
speed of light rays. It does not matter whether c+ is cal- 
culated with respect to the physical or the conformally 
rescaled metric because the null cone structure is invari- 
ant under conformal rescalings. We depict in Fig. 3] null 
curves on our grid. We see that the difference in the ini- 
tial time shift required for the synchronization is almost 
negligible for observers sufficiently far away. The reason 
for this is the fact that surfaces of a future hyperboloidal 
foliation follow closely outgoing null surfaces. Note that 
in a Cauchy-type foliation, the synchronizing shift in time 
would be calculated with respect to a far away observer 
at a finite distance. The required shift would then be on 
the order of the distance of that observer. 



t/m 




Fig. 4: Grid diagram depicting outgoing null curves on the 
grid. To synchronize observers from 100m to infinity, we need 
to take care of the shift in time given by the height of the thick 
curve. The diagram shows that the required shift is negligible 
in comparison to the time scales we are interested in. 

The next step for the construction of a physically moti- 
vated time parameter would be the reparametrization of 
the time coordinate so that the measurement by a given 
observer is plotted in proper time. The relation of the 
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proper time of a Killing observer at f = R to the coordi- 
nate time is given by 



As 



V '-gttdt 



2m 



1/2 



At. 



We see that at infinity proper time corresponds to coor- 
dinate time. The rescaling makes only a small difference 
for observers far away. 

A transformation of the results presented in this article 
according to the choice of such a physically motivated 
time coordinate does not change the conclusions for far 
away observers. Therefore, and because of simplicity, 
the plots will be given in the local time coordinate of 
the numerical solution. Note, however, that we do not 
plot the decay rates for observers near the horizon in our 
diagrams although our coordinates do allow this as well. 



4- A remark on finite distance extractions 

In numerical calculations where an artificial timelike 
outer boundary is introduced, one extracts gravitational 
radiation information from the numerical solution at fi- 
nite distances away from the source. When the extracted 
information does not vary in an unexpected manner with 
extraction radius or algorithm, one suggests that the cal- 
culation delivers a reliable answer. This practice can be 
regarded as safe as long as one has a theoretical under- 
standing of the expected radiation signal which delivers 
tests of the result, such as the peeling behavior, conserva- 
tion of mass or waveforms from approximation methods 
to compare with 0, [27j ■ A naive interpretation of finite 
distance extractions without theoretical knowledge, how- 
ever, can be misleading in terms of results as well as error 
estimates. 
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Fig. 5: Observer dependence of the tail decay rates at sub- 
sequent times. The thick curve corresponds to T — 20000m. 



The decay rates at late times provide an example for a 
case which demands caution. In Fig. Owe present a dif- 
ferent way of looking on the same data as given in Fig. [3] 
Here, decay rates at different finite distances have been 



plotted on subsequent time steps. The thick curve corre- 
sponds to the observer dependence of the decay rates at 
late times. We see that a direct Richardson extrapola- 
tion based on finite distance decay rates would not only 
deliver a wrong result, but also a wrong error estimate. 



5. Scalar field with angular momentum 

The difference in the expected tail behavior at null 
infinity and at finite distances becomes larger with 
At finite distances one expects the exponent to be 
p = —21 — 3, whereas at null infinity p = —1 — 2. In Fig. [6] 
we plot the decay rates for I = I after the quasinormal 
mode ringing is over and the tail part of the solution re- 
mains. In this case, there is no conjecture of an optimal 
pointwise decay estimate yet available in the literature, 
so we only plot the result of the numerical calculation. 



p(t,R) 
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Fig. 6: Decay rates for a solution to the scalar wave equation 
on a Schwarzschild background with I — 1 plotted for the 
same observers as for Fig. [3] This plot covers a shorter time 
than Fig. [3] due to the faster fall-off of the field. 

The qualitative behavior depicted in Fig.[5]is similar to 
the one depicted in Fig.[3]which suggests that an estimate 
similar to (| 1 3|) with modified exponents can be expected 
also for higher multipoles. 



B. Tails for the Yang-Mills equation 

The conformal method can be applied to study other 
fields as well. The basic requirements are that the equa- 
tions transform in a well-defined manner under conformal 
rescalings of the metric and the fields have a suitable 
asymptotic behavior. Yang-Mills equations are confor- 
mally invariant so that we can study them directly with 
the conformal method. 

The tail behavior for Yang-Mills equations have been 
studied in [f~lT | using linear perturbation analysis. They 
predict p — —5 at future timelike infinity and p = —3 
at future null infinity. Bizon et. al. showed in Q that 
the linear perturbation theory applied in [Tl| does not 
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capture the correct tail behavior for Yang-Mills equations 
where nonlinear terms appear (jTTJ) . Instead they argued 
and numerically demonstrated that the tail exponent at 
timelike infinity is p = —4 both for Minkowski and for 
Schwarzschild spacetimes. 

p(t,R) 
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Fig. 7: Decay rates for a solution to Yang-Mills equations 
on a Schwarzschild background. Each curve depicts the 
function p(t, R) calculated as in (|14|l for the following values 
R/m = (oo, 40000, 20000, 10000, 4000, 2000, 1000, 400, 200, 100). 



We plot in Fig. [7J the observer dependence of decay 
rates for a solution of the Yang-Mills equations. Near fu- 
ture timelike infinity a decay rate of p = —4 is observed 
in accordance with [5j]. One might expect that the pre- 
diction of at future null infinity will also be wrong 
by the same rate which indeed seems to be the case. As 
seen in Fig. [71 the numerical calculation suggests a decay 
rate of p = —2 for Yang-Mills fields at null infinity on a 
Schwarzschild background. 



V. CONCLUSIONS 

We investigated the asymptotic behavior of spherically 
symmetric solutions to scalar wave and Yang-Mills equa- 
tions on a Schwarzschild background. The investigation 
develops an idea presented in [401 ] where, among others, 
decay rates have been used to estimate the validity of 
the concept of null infinity for predicting detector signals. 
Our investigation focuses on the observer dependence of 
decay rates. In the case of scalar fields with vanishing 
angular momentum, a powerful mathematical estimate 
presented in [42|, [43[ could be compared with a numerical 
solution. For scalar fields with angular momentum, no 
such estimate is yet available, but the numerical results 
suggest that a similar estimate with modified exponents 
should be valid. For Yang-Mills fields the decay rates 
near timelike infinity and null infinity have been studied. 
The breakdown of linear perturbation theory for predict- 
ing decay rates near timelike infinity in case of Yang-Mills 
fields as discovered in [B| could be numerically supported 
also for decay rates along null infinity. 



On the example of the above studies, we emphasized 
the astrophysical relevance of null infinity. The exam- 
ples are intended to contribute to the clarification of the 
question of to what extent null infinity can be regarded 
as a good idealization for far away observers. The stud- 
ies also deliver a cautious remark against naive inter- 
pretations of finite distance extractions. One needs to 
be careful in wave extraction not only due to gauge is- 
sues in extraction methods and influences from the outer 
boundary, but also due to the fact that the extraction 
radius in standard numerical calculations is very small in 
comparison with astronomical distances. While the tail 
is certainly a special feature in that respect, there might 
be other features in the radiation signal that show a non- 
trivial dependence on the extraction surface. Advanced 
numerical studies including null infinity should be car- 
ried out to clarify this issue as has been initiated within 
the characteristic approach 24 1. 

The numerical studies of null infinity presented in this 
article have been performed by solving hyperboloidal ini- 
tial value problems 0, [2j|. This method seems to be 
quite powerful to study radiative properties of test fields 
at null infinity once a suitable class of gauges has been 
chosen 0, [T3, H3| • It also promises to be more general 
than the characteristic approach due to the large free- 
dom involved in the choice of hyperboloidal gauges. In 
comparison to the standard initial boundary value prob- 
lem, the numerical implementation is simpler because one 
does not need to mathematically construct or numerically 
implement boundary conditions and data. In addition, 
the conclusions that can be drawn from the numerical 
simulation are more powerful as radiative properties of 
fields can be studied all the way up to null infinity. 

In future work, the hyperboloidal initial value prob- 
lem should be further studied, for example, for different 
test fields or on different backgrounds such as on Kerr 
spacetime using a conformal fixing gauge as presented in 
[461 ] . The most important challenge in numerical stud- 
ies of hyperboloidal initial value problems is to devise a 
successful numerical implementation of the conformally 
compactified Einstein equations including null infinity. 
While there are some recent ideas on how to treat a hy- 
perboloidal initial value problem for the Einstein equa- 
tions directly [l], H3, [H, HU , they have not yet been im- 
plemented numerically for interesting cases. Studies of 
test fields might give further insight into this challenging 
problem. 
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